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The aim of this comment is to discuss and explain several aspects related to the derivation and application of the 
RPA theory in relation with a recent paper of Dinh Dang, there RPA theory has been applied to the two-level 
pairing model. We have some criticisms and remarks. We will start with the RPA in the boson formalism as used by 
N. Dinh Dang We notice that Dinh Dang uses a quite unconventional and non systematic boson expansion. Let 
", us restate the standard boson mapping in the case of the two level pairing model 

O: H=^{N2^N,)-gnY.A]A,,, j,f = l,2. (1) 
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To lowest order one have 



m>0 
V " ,„>0 



, where &] , bi are ideal Bose operators. The occupation number operators obey the exact relations 

o ■ 



^2=5Z4ma2™ ^ 2blb2. (3) 



Dinh Dang now uses the particle number condition (Ni) + {N2) — (N) = N = 2il, which holds if, in the absence of 
I interaction, the lowest level is filled. From the above number condition, one obtains 

§ : (^'I^'i) = {blb2). (4) 

Dinh Dang deduces from this relation that bi — b2 — b, b\ — it is assumed that it is approximately 

valid to replace the two ideal bosons 61 and 62 by the single one b. In the first part of this work we want to study the 
; ^ ■ validity of this single boson approximation. Keeping the two bosons the pairing Hamiltonian is given to lowest order 
_C3 ; by 

H - ~^2{n - b\bi) + ^2blb2 - gn b\hj 
= -e(r! - b\bi) + eblb2 - g^{blb2 + bib\ 

+ blb\+bib2). (5) 



In the single boson approximation we have 



H = - -2{n - b^f b) + -2b^b- gn{b^b + bb^ +b''b'' + bb) 
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-e{n - b^b) + eb^b - .gSl(fo^6 + bfet + fot^t + ;,^)_ 



(6) 



Both Hamiltonians can trivially be diagonalized with the help of an RPA (Bogoliubov transformation) among the 
bosons. We obtain for |5j 

H = Vlib\bi + Vl2b\b2 (7) 

with, 



^2 = 5 + VsTi^e + 5(1 - 2f^) (8) 

and for © 

H = iub^'b (9) 

with, 



^ = 2eWl-^. (10) 



In Fig.n]we have traces these different eigenvalues as a function of y = gil/2e. Also shown are the exact values for the 
chemical potentials = ±{Eq^^ — Eq) where Eq are the exact groundstate energies obtained by diagonalization 
of the original pairing Hamiltonian. The reason why we compare with is given by the fact that the eigenvalues 
of standard pp-RPA have to be identified with these chemical potentials (see e.g. [ilSBI)- We also see that the 
eigenvalues ^1,^2 follow, at least for small values of V, quite closely the exact values. On the other hand the single 
boson approximation yields a completely erroneous result which seems to have nothing to do with the exact solution. 




FIG. 1: Excitation energies Eexc in the non superfluid region as function ofV = Q/2e described in the text, for particle number 
N = 8. The spin of the levels is J = 7/2. The results refer to to exact calculations = ±(i5^*^ — Eq^) (double-dot dashed 
line, dotted line), standard pp-RPA Eqs.JHJ (solid line, dashed line) and pp-RPA in the single boson approximation Eg. llOII 
(dot-dashed line). 



Let us now comment in the superfluid phase. We note that in the superfluid phase we can calculate in standard 
QRPA the contribution of the q-term (see Eq. (7) in y|) of the Hamiltonian. Consequently, there is no reason to 
neglect this term as discussed by N. Dinh Dang The standard QRPA matrices are well-known and given by 



= 2(S, +2<7,y)%, +d,y (11) 
1 

n 



Bn' = 21- -&,r (12) 



where the quasiparticle energy is Ej — y (cj — gvj — fj) + and the gap A, as calculated from the BCS equation, 
includes the self-energy term. We have used the same notation as indicated in the work of Dinh Dang. We remind, 
shortly, that in this case, the gap is given by 



together with 



A^XIg^n^-^, (13) 



' 2 V 2gnJ 



(14) 
(15) 
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where ^ is defined as ^ = 2eQ,/{2Q. — 1). Using the latter relations we can write 



and we can calculate the different contribution of each term Ej, qjj' , djji, and hjjr 



Therefore, explicitly, the matrix elements are given by 

A2 A2 



2gn 

A2 A2 

2gn^ ^ 2^ 
A2 



and the positive eigenvalues of the RPA matrix are given by 



(16) 



A_.u,v,u,v,^ulvl = ^ (17) 



i?i=i?2=5f^+^;gn=g22 = -^ (19) 

A2 A2 
dii = = -g^l + Tr-F:'^di2 = ^21 = -ttt^ (20) 
Zgil Zgil 

^2 , , .gO A2 



All = A22 = gn-^^ + — (22) 



A12 = A2i = -^ (23) 



Bn = B22 = -TTZ?^ + T7Z?^ (24) 



r!i=0, 1]2 = 2A-^1-^ (26) 

in agreement with the result found by e.g. Hagino and Bertsch ^3j. In the case of the Fermion formalism of Dinh 
Dang, to obtain the RPA matrix elements, we have to begin with (see Eqs. (57-61) in [l|) 



Ajf = 2 {Ej + 3qjj>)Sjj, + djf 
1 



Bn' = 2 1--%' Uj/- (27) 



where the factor 3 is different from the correct factor appearing in H12|) . Exphcitly, in the case (|27() . the QRPA 
matrices are given, as follows 

= = ^29) 

= = - w + ^ ^''^ 

= 521=.^^-^ (31) 

where the gap A as calculated from the BCS equation includes the self-energy term. The RPA eigenvalues are given 
by (see Eq.(65) in Q) 



When we neglect the q-term (see Eq. (7) in we obtain 



^^2 = 2AJ(1-A) ix^-^\. (33) 



A,. = A,, = (35) 
A2 A2 

^- = ^-=^"-^ (37) 



what leads to different RPA eigenvalues of that given by Dinh Dang for this case. However, if we divide the q-term 
by a factor of 2, we obtain 

A2 

All - =gf^+^ (39) 
A2 

Ai2 = ^21=-— (40) 



2g?L 

A2 A2 
2^ ^ 2^ 



Bll = ^22 = -7^ + 777?^ (41) 



A2 

Si2 - S21 - .9^^ - (42) 



which produces the following RPA eigenvalues 



which coincide, exactly, with the solution given by Dinh Dang in the case where he neglects the q-term. We see that 
in this case the Goldstone theorem is violated and therefore the particle number symmetry is not restored. 

Dinh Dang also treats within the Fermion formalism several other approximations which do not produce the 
Goldstone mode at zero energy. One does not very well understand the aim of these considerations, since any way it 
is well known that the Goldstone mode should come at zero energy to restore the particle number symmetry 

which otherwise is violated. 

As a last point we would like to mention that Dinh Dang is superposing addition and removal modes in the non 
superfluid phase in Eq. (75). This superposition couples, simultaneously, N±2 states and therefore, obviously, violates 
particle number conservation already in the non superfluid phase. On the other hand it is well known y, |^ that 
pp(hh)-RPA in the non superfluid phase perfectly respects particle number symmetry : the addition mode involves 
amplitudes {N + 2\a^a^\N) and the removal mode {N — 2\aa\N). None of these amplitudes violates particle number. 

In short, in this comment we pointed out a number of shortcomings and inconsistencies in the work by Dinh Dang 
involving pp-RPA and QRPA in a solvable model. 
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